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Donor-Acceptor One Step Energy Transfer via Exchange 
Conpling on a Fractal Lattice 

C. L. Tang, M. A. El-Sayed* 

Department of Chemistry and Biochemistry, University of California 
at Los Angeles, CA 90024, U. S. A. 

ABSTRACT 

Temporal behavior of donor intensity, I D (t), resulting from one step 
donor-acceptor electronic excitation energy transfer process via exchange 
mechanism is calculated on fractal lattices, with discrete dilation symmetry, of 
Euclidean dimension, d • 2, and fractal dimension, D, ranging from 1.99 to 1.0. 
Ip(t) is fitted to the approximate equation of Klafter-Blumen (K-B) which is 
useful in fitting experimental results to determine fractal dimension from the 
slope of the expected straight line obtained by plotting In (- In Ip(t)) va. 
In (In** t). The result for fractal lattice with D ■ 1.99 indicates that the 
approximation is appropriate for the time range used in our calculations. The 
results for different fractal lattices also show that the K-B equation indeed 
gives a straight line for structure when D/d is not much smaller than unity. As 
this ratio decreases, deviation from the expected straight line results and an 
oscillatory behavior is observed. From the oscillatory characteristics, 
strnctnral information (i.e., fractal dimensionality, and geometrical parameters 
characterising the fractal lattices) as well as the molecular interaction 
parameter y, characterising the distance dependence of the exchange interaction. 
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INTRODUCTION: 

Fractals 1 , structures with dilation symmetry, have attracted a great deal 
of recent attention due to their nsefnlness in describing disordered systems. 
Processes occnring in systems snch as polymers 3 * 3 , proteins 3 '* and on surfaces 3 
have been discussed in terms of fractals, as are processes such as crystal 

£ 7 g 

growths dielectric breakdown* turbulence and cbaor • Fractals hare also been 

q 

used to describe the diffusion of liquids into porous media . 

Three different dimensions are at least required to define a fractal 1 ®* 11 . 

The first is the Euclidean dimension, d, in which the structure is embedded. 

The second is called the fractal dimension 1 , D. This describes the dependence 

of the number of sites N(R) on the distance R, through the relation (N(R) - R®). 

10 11 e 

The third dimenaion is the spectral or fracton dimension*'"*- 1 ’, d. which governs 

the random walk and relaxation processes and determines the density of states of 

the structure. The spectral dimension has been previously discussed in 

electron-spin relaxation studies in protein 3 ** and triplet-triplet annihilation 

studies in mixed molecular crystal 13 * 13 . 

Recently, studies of the one step electronic energy transfer have been 

discussed, both theoretically 1 * and experimentally 13 , in terms of the fractal 

dimension. Theoretically, Zlafter and Blumen 1 * (K-B) extended the equation 

ltf 19 

describing the time dependence of a donor intensity derived previously to 

fractal systems for one step energy transfer via dipole-dipole and exchange 
interaction. The derivations of these equations imply continuous dilation 
(i.e., the open fractal structure contains holes of any sixe). In our earlier 
work 11 , we have examined the fit of the calculated time dependent donor 
intensity for one step dipole-dipole trapping prooess on the generated fractal 
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lattice with discrete dilation symmetry. A structure possesses dilation 
symmetry if it remains equivalent to itself when one changes the unit length by 
a factor b. If b can be any number larger than one for thia structure, the 
dilation symmetry is continuous. Contrary, if b can only be some discrete 

values, the dilation symmetry will be called discrete. Structures with discrete 

dilation symmetry do not possess translation symmetry and contain holes of 

different sizes as classical fractals. Due to the discrete property, the law. 

N(R) ~ B®, which governs the varations of the mass N of the structure as 
function of the distance R, applies on the discrete system only if N(R) is 
averages over an interval from R to bR. For these reasons, these structures 
exhibit characteristic microenvironments which are repeated at different length 
scales. Therefore, calculated time dependent donor intensity plotted in the 
In -(In I(t)) versus (In t®^) showed deviation from the straight line predicted 
from the K-B equation 1 * which is based on the continnous dilation symmetry. Hie 
deviation showed an oscillating behavior with a constant periodicity. The 
amplitude of the oscillation increases as the ratio D/d decreases. Also, from 
the dilation symmetry criterion, the effect of changing the acceptor 
concentration is equivalent to a change of time scale for the measurement. This 

behavior is indeed observed in the plots with different acceptor concentrations, 

17 

P^, used in the calculation 

In this work, we will examine another kind of interaction, exchange 
interaction, which is the mechanism reponsible for triplet energy transfer in 
many aromatic compound systems . The effect of the variation in the 
microenvironment of the donors on the fractal structures with discrete dilation 
symmetry is checked for the exchange interaction. The temporal behavior of 
donor intensity undergoing a one step energy transfer via exchange interaction 
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on a fractal structure is calculated. Due to the difference in the system for 
which the K-B equation’®’* is derived (system of continuous dilation)* deviation 
froa straight line behavior is expected. Ve need to find answers to the 
following questions: a) can we observe the oscillatory behavior for 
exchange-type transfer siailar to that observed for dipolar tx_nsfer with the 
same parameters (i.e., fractal dimension D, acceptor concentration P^)? b) do 
the scaling properties of the donor emission intensity that relate to changes in 
acceptor concentration and of the tiae scale exist or not in exchange type 
transfer? c) for exchange interaction, is there a sealing property relating the 
interaction parameter, y, (which measure the range of the exchange interaction) 
and the time scale of the transfer? 

The calculation is carried out on several simple fractal structures all 
having Euclidean dimensionality d ■ 2 but with fractal dimensionality D “ 1.99, 
1.75, 1.5, 1.25, and 1.0. The acceptor concentrations end the interaction 
parameters y are varied to examine the dilation symmetry criterion. Like the 
results of dipole-dipole interaction, calculations show that the K-B equation is 
useful for most structures. A straight line relationship was obtained from 
which the correct D was determined for structures with D/d not largely different 
from unity. However, as D/d decreases, oscillatory behavior appears with 
characteristics relating to the different geometrical parameters describing the 
fractal structure. The dilation symmetry criterion is then examined and 
indicates that these properties may give us another approach to determine 
whether a system has the dilation symmetry property just by varying the acceptor 
concentration and/or changing the interaction parameter y by aelection of 
different pairs of donors and acceptors. However, due to the short time 
intervals and concentration usually employes in experimental studies, the 
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oscillatory behavior might be difficult to observe. 


THE MODEL: 


We consider only one step donor to acceptor energy transfer via exchange 
interaction. Donor-donor energy transfer is exclnded. Furthermore, energy 
levels of the donor and acceptor and the temperature are such that no acceptor 
to donor back transfer is allowed. 

Consider the case of an excited donor located at position Bp which is 
surrounded by a statistical distribution P A (B^, B^) of acceptors at position B* A . 
These acceptors can trap donor excitation with transfer rate V(H) which depends 
on the relative distance B « j B^-Rpjbetween the acceptor and donor. The decay 
curve of the donor excitation due to energy transfer follows the equation: 


I n (t) « exp 


-(J d* A P A <VK A > (1-exp ( -tW(B)))) 


where dB^ is the volume element over which the interaction is carried out. 
Equation (1) has been used previously assumming all donor positions are 
equivalent. For cases in which restricted volume exists, the acceptor 

•A afc 

distribution function P a (Rq«& a ) can be varied for different donor sites. For 
this reason, sveraging over the donor position is necessary in our lattice 
model. The averaging of equation (1) over the position of the donors with 
distribution P D (B^) leads to: 


I D (t) -J PjjCBp) dfi D exp -(JdB A P A (B^.B A )(l " exP -tW(B))) 
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The transfer rate f(R) from the donor to the acceptors for exchange interaction 

is 

W(R) - exp (-y(d-R)) (3) 

where y (interaction paraaeter) is a neasnre of the range of the exchange 
interaction, d is the nearest neighbor distance and R is the relative distance 
between the donor and the acceptor. Substituting Eq. 3 into Eq. 1, and 
integrating by parts, Eq. 3 can then be expressed as 2 ^ 

I D (t) - exp( j exp(-X) (In Z - In I) D dX) (4) 

where X *» exp ( y R), and Z “ (t/t). 

The term (In Z - In X) D can be expanded in power series for non-integer D, 

In X D(D-l) In X 

(In Z - In X) U - In Z (1 - D (-) +-(-) 2 + ...) (5) 

In Z 2 In Z 

end can be approximated as In® Z (the first tera) when the remaining terms are 
auch smaller. The approximation is valid only if 

In Z 

- > > 1 , ( 6 ) 

D (In X) 

which means the approximation is good when the experimental measurement is 
carried out at time much longer than the near neighbor transfer time. The donor 
intensity function can be simplified to 2 ** 
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In I(t) ~ ltP t (7) 

which will b* used to check our calculation reaulta. 

To get discrete dilation symmetry, we nae an idea similar to that used by 
Fournier in his galaxy model*. Details of generating the lattice is described 
in our earlier paper*^. Figure 1A shows the effect of the first iteration on 
the generation pattern. The dilation factor b is determined via the relation*^: 

b - N o */D (8) 

where it the number of points in the generating pattern, and equals to 4 in 
this calculation. Two examples with two different fractal dimension. D ■ 1.7 
and l.S. are shown in Fig. IB, and 1C respectively. One can see, from these two 
examples, that smaller the D value is, the more open the structure becomes 
(i.e., with larger holes). 

Due to the fact that a donor located near a large hole will have different 
environment than a donor located away from the hole, one has to average over all 
the different possible sites. However, in order to avoid the problem arising 
from the finite size effect of the fractal (due to the finite number of 
iterations used in the calculation), the donors were randomly distributed on the 
central area of the lattice, while the acceptors were randomly distributed on 
the whole fractal structures. The number of sites within the central area is 
sufficiently large to reveal the geometric characteristic of the discrete 
fractal structure. Also, as mentioned in the result section, the finite size 
effect has been examined by the dilation symmetry criterion: a) the calculated 
time dependent intensity I(t) with y ■ a can be mapped to the result with y * na 
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when changing the tine scale by the relation ahown in Eq. 9. b) the calculated 
I(t) with acceptor concentration, P^, can also be mapped to the result with 
P a /N c when increasing the transfer rate by exp(yb) according to Eq. 10. The 
number of iteraction we use to build the whole fractal is 5, which leads to a 
total number of acceptor sites of 4096. The number of iterations for the 
central area is 4 and the corresponding total number of donor points is 1024. 


RESULTS: 


The time dependent donor intensity function, Ip(t), is calculated by using 
Eq. (2) at different time t, for different acceptor concentrations, P^, for 
different fractal dimensions D, and for different interaction parameters, y, for 
exchange interaction. A DEC VAX 11/780 computer is used. Numerical integration 
is carried out on the generated lattice as described in the model section. 

The results of time dependent donor intensity function are plotted in 
accordance with the equation** assuming continuous dilation symmetry (see Figs. 
2 and 3). Deviations from the predicted equation due to the effect of 
'openness' in the fractal structure as the function of its dimensionality (i.e., 
the degree of compactness), as well as the effect of interation range are 
examined. Ve plot either logarithm of I(t) versus logarithm of t raised to 
power D or the logarithm of the negative logarithm of I(t) versus logarithm of 
In® t. Either plot should give a straight line as equation (7) predicts. 


A: THE D (FRACTAL DIMENSION) DEPENDENCE: 


One can see from Figs. 2 and 3 that deviation from the behavior predicted 
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foi the case of continuous dilation syaaetry occnra. As the fraotal diaension 
decreases, the deviation becoaes larger. These deviations are oscillations 
which increase with tiae in Fig. 2, and are of constant amplitude as shown in 
Fig. 3. The derivation of Eq. 7 has need the assumption of long tiae behavior 
to siaplify the equation. This approxiaation can be verified by the result of D 
■ 1.99 where the anplitude of oscillation should be very saall, and the 
deviation froa the linearity should coat froa the long tiae epproxiaation. In 
Fig. 2, straight line behavior is observed, which suggests that use of the long 
tiae approximation in the calculation is valid. 

B: THE EFFECT OF y (INTERACTION PARAMETER) ON TIE TIME DEPENDENCE Of THE DONOR 
INTENSITY: 

The result of the donor intensity decay function for different interaction 
parameter y for exchange interaction (i.e., D - 1.0, and P^ - 0.05) is shown in 
Fig. 4. One can notice the differences between different y's are the scales of 
the intensity and time. Froa the property of dilation syaaetry, the calculated 
intensity of y ■ a can be mapped to y ■ u by the transformation: 

I D (t.y - a) - I D (t 1/n .y - na). (9) 

This relation can be used to check the finite size effect of lattice sites 
(minimum length scale for the finite nearest neighbor distance and maximum 
length scale for the total number of sites in the lattice, 40 96) due to tiae 


limitation for the computer calculation. If the number of iterations does not 
purturb the result, one should expect no change in the time behavior if 
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nultiplying the value of y by n, and decreasing the tiae scale by the power of 
1/n. This is observed in Fig. 4 (cnrves A and B), where the value of y ia 
changed froa 5 to 10 and it is equivalent to change the In** t scale to 1/2 of 
ln D t (i.e., 40 to 20 for D -1.0). 

C: ACCEPTOR CONCENTRATION AND THE TIME DEPENDENCE OF THE DONOR INTENSITY: 

In parts A and B of Figs. 2 and 3, the concentration of acceptor, P^, is 
0.03 and 0.012S, respectively. One can see the difference between thea ia 
aainly the scale of the intensity. The shape and the curvatnre of these two 
plots are exactly the saae (i.e., one can superiapose one plot on the other). 

Like Eq. 9 in part B of this section, there is siailar dilation syaaetry 
criterion for the relation between tiae and concentration, 

I D (t * P A* ” I D (t*exp(yb). P A /N 0 ). (10) 

where b is the dilation ratio (Eq. 8), and N q is the number of points in the 
generating pattern. If the number of iterations does not purturb the result, 
one should expect that by dividing the acceptor concentration by N q (N q is equal 
to 4 in this generating pattern) and increasing the transfer rate by exp(yb), to 
get the same tiae behavior. This is due to the fact that increasing the 
concentration of the acceptor by a factor 4 (N q ) *nd keeping the tiae scale 
constant is equivalent to holding the concentration constant (i.e., the sane as 
decreasing the length scale by a factor b), and decreasing the tiae scale by a 
factor of exp(yb). This infers that Figs. 2B and 3B are the short tiae part of 
Fig. 2A and 3A. 















PAGE 12 


DISCUSSION: 

la oar earlier paper*^, ve have shown that fractal latticea with discrete 
dilatioa symmetry create aaother kind of exclnded volume (i.e., holea) which ia 
characterized by differeat shapes of geaerating pattera of the lattice. Vhea 
these holea hecoate too important to be neglected, deviations from the temporal 
behavior predicted for fractals with contianons dilation symmetry are observed. 
In earlier calculations*^ for one step trapping procesa with dipole-dipole 
interaction, we have observed that deviation from the equation** ahowed 
oscillatory behavior with constant periodicity in the In t ® !!> scale. The 
amplitude of the oacillation dependa on the ratio of D/d. The amaller the ratio 
ia, the larger the amplitude wonld be. 

In thia work, the effect of interaction type is examined. Since the 
difference between the dipole-dipole and exchange interactions is the range of 
transfer process at a given time (i.e., the time dependent function for the 
transfer rate), we would expect to see the similar oscilletory behavior with a 
constant periodicity in the In In® t scale deviate from the predicted equation 
derived for continuous dilation symmetry like the dipole-dipole trapping ease. 
As shown in Fig. 3, the oscillation is indeed observed in the calculation. The 
reason for this oscillatory behavior is the same as that mentioned for the 
dipol-dipole case,*"* except the x-scale in Fig. 3 has been replaced from In t®^® 
to In In® t because of different time dependent function form of the transfer 
rate for each interaction type. The effect of changing the interaction 
parameter y in the exchange interaction is examined. In Eq. 3, we can see that 
for the same transfer rate T, as the value of y increases, the value of R 
<R * | Sp - R^| ) decreases. Therefore, the larger y is, the shorter the 
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distance the transfer probes. Because the total number of acceptors within the 
distance B is proportional to B® (the relation is exact only when averaged from 
B to bB as shown in the introduction section), the donor intensity would show 
smaller decay as y increases. This behavior is shown in Fig. 4. 

As mentioned in the result's section* we have used the dilation symmetry 
criterion to check our calculation for the finite size effect. The dilation 
symmetry criterion relates the change of the acceptor concentration to the 
change of the time scale via Eq. 10. For exchange interaction, there is another 
criterion that correlates the change of the interaction parameter y to the 
change of time scale via Eq. 9. In principle, this means that a knowledge of 
the time dependence of the donor decay at a given acceptor concentration (or 
interaction parameter) should enable us to predit the time dependence of another 
acceptor concentration (or interaction parameter) and/or time scale. Both 
correlations can be observed in Figs. 2 and 4. where the acceptor concentrations 
are 0.03 and 0.0123 (for Fig. 2) and the range parameters are 1., 2.. 3.. and 
10. (for Fig. 4). One can see the only difference between them is the scale of 
the intensity and the time. The shape and the curvature of these plots are 
exactly the same (i.e.. one can superimpose one plot on the other). In other 
words, the change of the acceptor concentration reveals that Fig. 2A is the 
short time behavior of Fig. 2B when one changes the time scale by a factor 
exp(yb), or Fig. 2A is the magnificstion of the upper left-corner of Fig. 2B. 
Similar results are observed for the chsnge in intersction parameter, curve A in 
Fig. 4 (y - 10) is the short time behavior of curve B (y ■ 5) when one chsnge 
the time scale t to t*^ (or 1/2 in the In® t scale, D - 1 in Fig. 4). 

All the properties mentioned above are results of the dilation symmetry of 
the system under study. Also, information about the structure of the lattice 


' */ - * O - V.* • - - • O •/ * ■ *' • 0 
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can be obtained fro* the oscillatory behavior. For example, the distanoe 
between two aaxiaa in the oscillation in Fig. 3A and 3B is 1.39 in the 
coordinate of In In** t for D equal to 1.0). This coordinate is eqnal to the 
logarithm of the length scale with power D (i.e.. In 1®). Therefore, the 
corresponding value of R between two maxima is eqnal to 4. which is the valne of 
b (i.e.. in this case. D - 1.). The change in the logaritha of negative 
logaritha of the intensity corresponding to this period is eqnal to In N q , where 
N o is 4 for the generating pattern we choose. Fntheraore. D still can be 
deterained froa the slope of the best straight line fit through the oscillating 
function. Finally, the dilation symmetry criterion can give ns another way to 
exaaine whether the systea investigated possesses dilation symmetry, 
experimentally, by chsnging the acceptor concentration and/or interaction 
paraaeter for exchange interaction by selecting different donor and acceptor 
pairs. 

In principle, the oscillatory behavior of the deviation from the predicted 
equation should give us a lot of information about the structure of the systea. 
However, in real systems, due to the limitation of the lifetime and the transfer 
time of the donor excited state, experimental determination of real behavior is 
difficult. For example, in triplet-triplet energy transfer between aromatic 
hydrocarbons (which occurs via the exchange mechanism ), the radiative lifetiae 
for the triplet state of these molecules is on the order of tens of seconds. If 
the transfer rate is too fast, experimentally, it would be difficult to follow 
the temporal behavior of the energy transfer process. Therefore, only a small 
tiae range (i.e., 10”* to few seconds) can be studied experimentally. Assuaing 

that the near distance is d and the transfer tiae x within this distance d is 

—12 —4 
~ 10 sec, then, the distance probed at 10 sec is about 4.7 d for y - 5 (as 


! 
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was inferred for aromatic molecule* in crystals ), and the distance probed at 
10 secs is 7.0 d. Crude estimation shows that the distance range probed in the 
triplet-triplet transfer process for moat aromatic hydrocarbons wonld be only 
about 2.3 d, and it wonld be very difficult to tell the real structure of the 
system studied from a simple one step energy transfer process. 

The existence of oscillatory behavior indicates that for some types of 
lattices which follow the equation: N(E) ~ tP, for the fractal structures, it 
is not necessary that the temporal behavior for trapping process do follow the 
equation derived from the continuous dilation symmetry. However, the equation 
would describe the temporal behavior well if the systems considered really 
possesses the continuous dilation symmetry and the oscillatory behavior would 
disappear. 

CONCLUSIONS: 

From the calculation of the temporal behavior of the donor intensity on a 
fractal lattice with discrete dilation symmetry, the K-B equation is found to 
describe most of the time dependent behavior for structures with D/d ratios 
close to unity. In these cases, a plot of In of negative logarithm intensity 
versus In of the logarithm of the time shows a straight line, and the resulting 
slope gives the value of the fractal dimension D. However, as the D/d ratios 
decrease, deviations from linearity become large, and oscillatory behavior 
observed. Fortunately, from the oscillatory behavior, information regarding the 
geometry of the fractal structure under study can be obtained. 

In principle, the results obtained for the cases of small D/d, for 
different acceptor concentrations and/or for different interaction parameters 
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should enable experimentalists to test whether the anomalous decay is dne to n 
real fractal space or to a regular structure with excluded volume. For fractal 
structure, it is possible, at least in principle, to characterise the fractal 
dilation syasetry experimentally by determining the value of and b from the 
oscillatory, temporal behavior of donor intenaity. Unfortunately, the ranges of 
time and concentration which can be studied experimentally in energy tranafer 
processes do not allow observing fnll oscillation. Instead, a part of the 
oscillation could only introduce errors in equation (7) when used to determine 
D. Furthermore, the value of D determined conld depend on the time interval 
plotted (i.e.. on different portion of the oscillation used). Urns, due to the 
experimental limitation imposed by the lifetime of the donor excited state, and 
the short range of the transfer of the exchange interaction mechanism, the 
determination of the real fractal distension from one step exchange energy 
transfer might be met with uncertainty. 
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Figure captions: 


Fig. 1: 

(A) Plot of the generating function of the fractal lattice (N q - 4) used in this 
study. The value of b for different D's are calculated by Eq. 9. (B) and (c) 
are the two typical fractal lattice structures obtained with this generating 
function for D*1.7 (scaling factor b “ 2.26) (B), and l.S (scaling factor 
b - 2.52) (C). 


Fig. 2: 

Plot of the calculated In of the donor intensity versus ln D t for different 
fractal dimension D. and for different acceptor concentration P^.— — — D*1 .99, 

-D»1.75, —-D-l.5.-D»1.25, --—D*»1.0. In part A, and B, the 

acceptor concentration is 0.05 and 0.0125. respectively, t is in arbitrary 
unit. The characteristic transfer time has been kept constant in Fig. A and B, 
but are different for each concentration. 


Fig. 3: 


Plot of the ln(-(ln(intensity)) versus In ItP (time) 
dimension D and for different acceptor concentrations. 


for different fractal 
The slope of the curves 
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should be equal to 1 in these coordinates if the fractals have continnons 
dilation symmetry assumed by ref. 14. The slopes calculated are oscillating 

with tine.-D- 1.99. -D-1.75,-D-l .5.-D-1.25,-D-l .0. 

(A) and (B) correspond to acceptor concentrations of P^-0.05, and 0.0125 
respectively, t is in arbitrary nnit. The best straight lines fit through all 
are parallel to each other. 


Fig. 4: 


Plot of the In (intensity) versus ln D (time) for different interaction 
parameters y. A: y-10., B: y“5.» C: y m 2., D: y m l . All the data calculated are 
with the same fractal dimension D (D-1.0) and acceptor concentration (P A «0.05). 
t (time) is in arbitrary unit, and is the same for all four curves. The curve A 
can be mapped to curve B by changing the time scale via Eq. 9 due to the 
dilation symmetry property of the generated lattice. The same procedure can 
also be applied to curve B, C, and D. 
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